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Abstract 

We perform extensive monomer-resolved computer simulations of suitably-designed amphiphilic 
dendritic macromolecules over a broad range of densities, proving the existence and stability of 
cluster crystals formed in these systems, as predicted previously on the basis of effective pair 
potentials [B. M. Mladek et ai, Phys. Rev. Lett. 96, 045701 (2006)]. Key properties of these 
crystals, such as the adjustment of their site occupancy with density and the possibility to heal 
defects by dendrimer migration, are confirmed on the monomer-resolved picture. At the same 
time, important differences from the predictions of the pair potential picture, stemming from steric 
crowding, arise as well and they place an upper limit in the density for which such crystals can 
exist. 



PACS numbers: 82.70.-y, 82.30.Nr, 64.70.D-, 83.10.Tv 



Soft cluster crystals (also termed bubble solids), i.e., solids where each lattice site is oc- 
cupied by several, interpenetrating particles, arise in systems governed by generic, purely 
repulsive interactions which remain finite even if particles fully overlap P-[3] . The completely 
different character of these novel cluster crystals, as opposed to hard cluster phases ^ [5] , is 
testified by the fundamentally distinct nature of the interactions between their constituent 
particles. The clusters are not stabilized by a combination of short-range attractions and 
long-range repulsions, which result in well-defined aggregates in the liquid that disappear at 
high densities [U E]. Rather, clustering is here a self-enhancing, cooperative phenomenon, 
where full overlaps with few particles are favored over costly partial overlaps with many 
neighbors. Thus, clusters remain stable due to their mutual repulsions with the surround- 
ing ones |6j. Cluster crystals are thus high-density phases, in which fuzzy agglomerates, 
already formed in the liquid, freeze into well-defined and spatially ordered groups at higher 
concentrations [3]. Predicted for certain macromolecular [Tl|2l[8] or atomic [T0Hl3] systems, 
such clusters form exotic crystals that show mass transport [HI [T3], unusual reaction to 
compression [31 EJ [12] and shear [T7] , and rich phase behavior with a cascade of isostructural 
transitions [TBI [19]. Further, they possess a unique rheological behavior, showing thixotropy 
and flow quantization [T7] . 

All hitherto performed studies of the properties of soft cluster solids have been based 
on the assumption of ultrasoft, pairwise additive, density-independent inter-particle inter- 
actions. According to theory, the necessary requirement for a system to display clustering 
behavior is that its pair interaction shows negative Fourier components [21 [20]. The pair 
potential picture is realistic when understood as an effective representation of certain types 
of polymeric macromolecules in not too concentrated solutions [2T1423] . However, the rele- 
vance of such an approach for experiments is here put into question since the formation of 
clusters is a high density phenomenon. Under such conditions, many-body effects, confor- 
mational deformations and crowding are expected to play a crucial role, which could either 
hinder or favor cluster formation. For ring polymers, the zero-density pair interaction has 
negative Fourier components; however, due to the shrinkage of the rings at finite densities, 
the clustering ability is lost [BJ. In this Letter, we show on a microscopic basis that soft 
cluster crystals can indeed be stabilized in a broad density regime for suitably designed 
macromolecular systems, motivating the experimental realization of such systems. 

Amphiphilic dendrimers constitute such a system predicted to show clustering [H [2]. 
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Here, we consider amphiphilic dendrimers of second generation with two central monomers. 
The dendrimers' 14 monomers are divided into two classes: the eight outermost monomers 
form a solvophilic shell which surrounds the solphophobic core build from the interior genera- 
tions' monomers. At sufficiently close separations, all monomers strongly repel each other via 
Morse potentials, while bonds between monomers are modeled by FENE springs [24J. In pre- 
vious studies, the zero-density effective pair interaction between various realizations of such 
dendrimers has been determined [1], showing the desired negative Fourier components |25] . 
Indeed, we then found the onset of clustering in the fluid phase of microscopically-resolved 
amphiphilic dendrimers [2] . Here we demonstrate the stability of the cluster solid phase (see 
scheme in Fig. SI of the Supplemental Material). Let p = N/V be the density of a system of 
N dendrimers in the volume V and Rg be the p = radius of gyration of such a dendrimer. 
An estimate, P/, of the freezing density p/ can be obtained on the basis of the infinite di- 
lution pair potential $eff('") [20] • Further, theory predicts the number of dendrimers per 
cluster, i.e., the occupancy iVocc, to be proportional to the density: iVocc oc p [2, 3J. Con- 
sequently, the lattice constant d oc {Nqcc/ pY^^ is predicted to be density- independent. For 
the dendrimers considered here, lattice site occupancies of Nocc = 3.5 (bcc) and Nocc = 3.2 
(fee) are predicted at the freezing density pjRg = 0.281. The lattice constants are predicted 
as a = 2.9Rg (bcc) and d = 3.6Rg (fee), leading in both cases to nearest-neighbor distances 
d = 2.hRg. 

These theoretical predictions are put to the test by performing NVT-Monte Carlo simu- 
lations. First, we delimit the range of values of A^^occ that can stabilize crystals: although in 
the effective description A^^occ can be arbitrarily large, in the monomer-resolved description 
local steric crowding is anticipated to put an upper limit to this quantity. To this end, we 
simulate disordered systems of dendrimers in a density regime around p/. We study the 
cluster size probability distribution P{Nocc), employing a cluster algorithm developed pre- 
viously [2]. Representative results are shown in Fig. [ijfor three different densities. Fig. [T]^a) 
shows P{Nocc) of a system at low density, pRg = 0.046, i.e., well below pf. Some small clus- 
ters are present but single dendrimers dominate: this system is a fluid of single dendrimers. 
Increasing the density to pRg = 0.138, Fig. [l|b), while still staying below p/, we flnd that a 
distinctive distribution of larger clusters has also developed. The sizes of these clusters are 
not uniformly distributed; rather a single cluster size, N^^fy is favored, in full agreement 
with theoretical predictions [20]. Here, the system forms a cluster fluid, as already studied 
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FIG. 1. Color online. The probability distribution P{Nocc) to find cluster size A'occ in (a) the 
fiuid state at pRg = 0.046; (b) the cluster fluid state at pRg = 0.138; and (c) the solid state 
at pRg = 0.323. The insets show simulation snapshots, where clustering can be seen from the 
crowding of the dendrimers' core monomers [black and blue (light gray)]. Shell monomers are 
shown in red (dark grey). Monomers are not drawn to scale. 

in detail in previous work, albeit for a different dendrimer model only showing clustering at 
higher occupancies than the present model [2l |26] . Finally, the density was increased above 
Pf to the value pR^ = 0.323, Fig. [l|c). Now, single dendrimers have almost completely 
vanished from the solution and have been replaced by clusters of a preferred occupation of 
N^cf = 4 — 5. Visual inspection reveals that this system has spontaenously developed local 
crystalline order, albeit forming a distorted crystal. This finding is corroborated by an anal- 
ysis of the radial distribution function of the clusters' centers of mass, which has developed 
a pronounced first peak centered at a value close to the predicted a and clearly separated 
from the second peak. System sizes, however, do not allow for a positive determination of a 
particular structure. Similar results were found for several other densities above pf, proving 
that the cluster liquid is unstable in this density regime. 

Spontaneous crystallization, however, does not provide information which crystalline 
structure is the one of lowest free energy; crystallites formed could exhibit the structure 
with the lowest barrier separating it from the fluid state ^71 12H] • We therefore artificially 
prepared perfect bcc and fee cluster crystals at various densities p, by placing Nocc ^ ^ocfip) 
dendrimers per perfect lattice site and then delimited their range of mechanical stability by 
performing NVT- and NPT-Monte Carlo simulations. Non-stable crystals melted within a 
few simulation sweeps into diffusive fluids, as manifested by a steady increase of the den- 
drimers' mean-square displacement and an accompanying loss of the ordered structure. The 
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FIG. 2. Color online. Simulation snapshot of a fee dendrimer cluster crystal of occupancy A'occ = 10 
at pRg = 0.438. (a) A full-monomer representation of the crystal. Clustering can be seen from 
the crowding of core monomers [black and blue (light gray)], (b) The dendrimers' centers of mass 
(small spheres) are found to 99% in the big, gray spheres around the cluster's centers of mass. 

mechanically stable systems, however, remained in the crystalline, clustered arrangement 
throughout the whole simulation runs. Fig. shows a typical simulation snaphot of a 
cluster crystal in the full-monomer representation of the dendrimers, where the clustering 
can be seen from the crowding of core monomers [black, blue (light gray)], which are sur- 
rounded by a cloud of shell monomers [red (dark gray); also cf. Fig. SI]. The clustering and 
the self-assembly of the clusters onto the lattice sites of a (fee) crystal become even more 
evident when plotting just the dendrimers' centers of mass. Fig. ^h), denoted there by small 
spheres. The irregular shapes of the various clusters are a manifestation of the instantaneous 
fluctuations of each dendrimer around its lattice site and show that the centers of mass of the 
individual dendrimers are not tightly bound to lie on top of each other. These findings are 
in full agreement with the theoretical picture, which predicts these oscillations to give rise 
to intra-cluster breathing modes that manifest themselves as optical branches in the phonon 
spectrum of the crystal jl8|. We quantify these fluctuations by plotting gray spheres, cen- 
tered around the clusters' centers of mass, and which contain the latter to a probability of 
99%. Spheres from different lattice sites do not overlap, confirming the prediction that the 
clusters are well-localized around their lattice sites |20j . 

The results of the simulations delimiting the range of mechanical stability are summarized 
in Fig. [3} revealing that a large, insular part of density-occupancy space allows for the 
formation of mechanically stable cluster solids. The minimum density to provide mechanical 
stability of the cluster crystals is found to be pRg = 0.225, independently of the occupation 
number. While the determination of the precise location of the freezing density p/ is not 
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FIG. 3. Color online. The mechanical stability of dendrimer cluster crystals for (a) bcc and (b) fee 
as function of density pi?^ and occupation number A'^occ- State points for which crystals remained 
stable are denoted by circles (•), those where crystals melted as diamonds (♦). The predicted 
freezing density is shown as vertical, dotted line at yoi?^ = 0.281. For high occupations and 
densities, systems transformed (a) or formed percolated networks (■) (see text). Systems for 
which the free energy was measured are circled in black. The dashed line plots the theoretical 
prediction Nocc{p), while the equilibrium occupation Nocc{p) (see text) is plotted as solid line. 



the subject of this paper, these results demonstrate that pf ~ pf. 

As can be seen in Fig. |3| the minimum occupation number leading to mechanically stable 
solids increases linearly with density for both bcc and fee. Systems with the lowest stable 
occupation, Nocc — 4 showed single events of lattice site hopping and subsequent healing of 
the resulting defect. However, due to the limited time-scales accessible in our simulations 
and the crowded systems studied, complex structural rearrangements and hopping events 
leading to long-time diffusion |TH [15] have not been observed frequently. Crystals of high 
occupations, A^occ — 14, further remained mechanically stable upon artificially introducing 
lattice defects like single clusters of lower occupancy, or even lattice vacancies or interstitials. 
Also mixed occupations, where not every lattice site was occupied by the same amount 
of dendrimers, remained stable, provided the average value of A^occ was chosen from the 
island of stability (see Fig. [3]), in full agreement with the theoretical prediction of cluster 
polydispersity in the crystal [T^ 120]. 

Crystals with even higher occupations, i.e., Nocc ^ 16 were not found to be stable: they 



transformed to distorted crystal structures of lower occupations by splitting clusters and 
halving the lattice constant, thereby taking these crystals back to the island of stability 
(Fig. [sj triangles). This shows the propensity of over-occupied crystals to adjust their 
occupancy towards an optimal value, a feature which will be confirmed quantitatively by 
free energy calculations in what follows. At high densities and rather low occupations, the 
lattice spacing accordingly shrinks and enters the length regime of the bonds between the 
dendrimers' monomers. Here, the crystalline arrangement of well-defined clusters also cannot 
be sustained: dendrimers stretch and are shared between two clusters, leading to disordered 
percolated networks of dendrimers (Fig. [sj squares). This phenomenon reflects that at such 
high densities the effective pair potential picture, which predicts an unlimited growth of 
cluster occupancy with density, A^occ oc p, breaks down qualitatively since steric crowding 
between the monomers prevents further cluster growth upon sufficient compression. 

Having established the existence of a broad region of their mechanical stability, it is 
pertinent to assess the crystals' thermodynamic stability and their optimal occupation at 
given density. To this end, we calculated the excess free energy per particle /3/ex = PF^^/N 
for the mechanically stable crystals with the help of thermodynamic integration as outlined 
in Refs. [HI [IS] ■ Here, the reference state is chosen as non-interacting dendrimers confined to 
the desired crystal structure by potential wells. Being the same for all crystal structures, the 
reference free energy (3Fid/N may be disregarded. Measuring /3/ex for several values of Nocc 
at given crystal structure and fixed density, we determine the equilibrium occupancy N^^^ 
as the one that minimizes the free energy, see Fig. |4](a) [19] . In qualitative agreement with 
theoretical predictions, N^^c increases linearly with density. Fig. |4]^b). However, there exist 
quantitative deviations between the two, manifesting the role of the many-body interactions. 
Whereas the pair-potential description predicts Nocc P, here there is a non-zero offset d, 
Nocc = cp + d. At freezing, the optimal occupation number from the monomer-resolved 
simulations turns out to be higher than theoretically predicted. This indicates that the many- 
body interactions, whose importance grows with occupancy, have the effect of strengthening 
the overall attractions between the soft particles at close approach. Since these properties 
are model-specific, we do not attempt a more detailed analysis here. However, we point out 
that the theoretical prediction of Nocdp) is found to be an excellent estimate for the lower 
limit in occupancy of mechanically stable crystals (see Fig. |3]). Indeed, these occupancies 
are sufficiently low for pair potentials to still be approximately valid, their degree of validity 
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FIG. 4. Color online. All data is plotted for fee (o) and bee (□) lattiees. (a) The free energy per 
dendrimer, /3/ex, as a funetion of A'occ at different densities, from bottom to top: pRg = 0.22, 0.28, 
0.36 and 0.44, showing for each density a different optimal occupancy, Kcc, that minimizes it. (b) 
The equilibrium occupancy, A'^^cc, as function of density pRg. (c) The equilibrium lattice constant 
a/Rg shrinks with increasing density, (d) The equilibrium free energy per volume, /SFc^Rg/V as a 
function of the density p, evaluated at the equilibrium occupancy Nocc{p)- Fee is found to be the 
thermodynamically stable crystal. 

diminishing as Nocc grows. In Fig. ^d) we show the free energy per volume of the equilibrium 
crystals as a function of density. We find that the fee crystals have a slightly lower free energy 
than bee for all densities considered. 

In contrast to amphiphilic dendrimers, ring polymers cannot form cluster crystals due to 
their drastic shrinkage upon compression [8] . It is therefore pertinent to study the conforma- 
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tions of dendrimers in the clusters of the thermodynamically stable crystals. The radius of 
gyration of the dendrimers shows little variation with density, due to the necessary balance 
between the exterior and interior osmotic pressure of clusters and the constraints imposed 
by the regularly branched dendritic architecture. Nevertheless, the configurations of the 
dendrimers change considerably upon increasing density, a fact manifested by a density- 
dependence of the most probable distance of shell monomers from their dendrimer's center 
of mass. While shell monomers can backfold to the dendrimer's core at low densities, the 
cluster formation at finite densities increasingly forces the shell monomers to move to the 
rim of the molecules |29j . 

Altogether, the cluster crystals discovered in our study are characterized by a hierarchy 
of self-organization at various interconnected levels. First, the individual dendrimers adjust 
their conformations to the thermodynamic conditions; second, the individual clusters adjust 
their population and are thus robust against cluster polydispersity and crystal defects; and 
finally at the macroscopic level, the flexibility of the cluster crystals is manifested through 
the mechanism of restructuring of non-optimally occupied cluster crystals, the hopping of 
dendrimers among crystal sites, and the adjustment of the associated lattice constant. Dis- 
crepancies between the fully microscopic and the effective potential pictures exist due to the 
fact that the strong steric repulsions of the monomers prohibit the ad lihidum growth of the 
clusters' population. Still, the salient features of soft cluster crystals have been confirmed. 
The range of stability of these structures is sufficiently broad to render them into materials 
of practical as well as fundamental importance. 

In summary, we have shown at a microscopic level that cluster crystals exist in suitably 
synthesized soft matter systems. More detailed computational investigations in both the 
effective and the microscopical description will be needed to uncloak the whole beauty of 
these intriguing systems and understand, for example, the nucleation of soft cluster crys- 
tals or the transformations of dendrimer crystals at high occupations and densities; this 
will be the topic of future investigations. Our results convey the clear message that an ex- 
perimental realization of cluster crystals employing simple and easy to synthesize dendritic 
macromolecules is possible. We thereby provide guidelines to achieve this goal, which will 
lead to the production of a new class of materials with most intriguing physical properties. 
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Pauh Institute. B.M.M. acknowledges the MFPL VIPS program (funded by BMWF and the 
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Supplemental Material for paper: 
Microscopically-resolved simulations prove the existence of soft 

cluster crystals 

Dominic A. Lenz, Ronald Blaak, Christos N. Likos, and Bianca M. Mladek 



THE DENDRIMER MODEL 



A generic view of the amphiphilic dendrimer architecture, of the effective pair potential 
and of the hierarchical assembly of the molecules into a cluster crystal is shown in Fig. 



SI We employ the dendrimer model already introduced before in Refs. [U [2] and based 
on Ref. [3]. We study amphiphilic dendrimers of second generation which have two central 
monomers. The dendrimers' 14 monomers are divided into two classes: the eight monomers 
of the outermost generation {g = 2) form the solvophilic shell of the dendrimer (index S), 
and all monomers of the interior generations {g = and g = 1) form the solvophobic core 
(index C), resulting in amphiphilic dendrimers. Interactions between any two monomers, 
separated by a distance r are modeled by a Morse potential p3], which is given by 

/3$;fr'(0 = { [e-"-('--'^-) - 1]' - l} ' /^^ = CC, CS, SS (1) 

where cj^i, denotes the effective diameter between two monomers of species fi and u. All 
energy scales employed are measured in units of (3 = {k-QT)'^, with T being the temper- 
ature and ks Boltzmann's constant. The Morse potential is characterized by a repulsive 
short-range behavior and an attractive part at long distances, whose depth and range are 
parametrized via e^j, and a^^, respectively. The core monomer diameter ac is chosen to be 
the unit of length for the model parameters given in Tab. [Tj 

The bonds between adjacent monomers of relative distance r are modeled by the finitely 
extensible nonlinear elastic (FENE) potential, given by 

2" 



R,, 



jju = cc, CS (2) 



where the spring constant K^j^j, restricts the monomer separation to be within a distance 
from the equilibrium bond length lfj_,y. 

We simulated more than 100 different realizations of dendrimers, determining the zero- 
density pair interactions via the algorithm outlined in [1]. We then chose the dendrimer 
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FIG. SI. Schematic representation of amphiphilc dendrimers of second generation, consisting of 
a solvophobic core (black and blue inner monomers) and a solvophilic shell (red outer monomers; 
monomers not drawn to scale). Several, interpenetrating dendrimers form clusters, which is ex- 
pected to be promoted by the dendrimer-dendrimer interaction, <I>eff(r) (bottom left), which favors 
clustering. Clusters can then freeze into crystals. Here, /3 = {kBT)~^ and Rg is the infinite dilution 
radius of gyration. 



with the deepest minimum in Fourier transform, which promises a high abihty of clustering. 
Its parameters are given in Tab. [T] 

We stress that we also carried out simulations of cluster crystals for other realizations 
of dendrimers, whose effective interactions are far less promising regarding their clustering 
behavior. While these systems are more difficult to access computationally due to their 
higher necessary equilibrium occupancies, the obtained results are in full accordance with 
the findings presented in the main part of the present contribution. 
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TABLE I. Potential parameters of the dendrimer considered in this study. The radius of gyration 
of an isolated dendrimer is Rg = 3.589(Jc. 



DENDRIMER STRUCTURE IN THE CLUSTER CRYSTALS 

In the main text, we refer to the adjustment of dendrimer sizes and conformations upon 
changes of the density p and the cluster occupation A^occ- Here we show quantitative ev- 
idence of this effect. The dependence of the dendrimers' radii of gyration on p and A^occ, 
Rg{p;Nocc), for the mechanically stable crystals in (cf. Fig. 3 of main text) is shown in 



Fig. S2 For reference, the zero-density value Rg is plotted as well. The equilibrium values 
of the radii of gyration R^lp), corresponding to the optimal occupancy {p,N^2c) according 
to Fig. 4 are displayed by the open symbols, revealing that their values are nearly density- 
independent and, interestingly enough, very close to the zero density value. The dependence 
of the dendrimer size on Nocc for a fixed density p is a monotonically increasing function, a 
feature that can be understood through the influence of two interconnected mechanism. On 
the one hand, as A^^occ grows, the dendrimers have to stretch to accommodate the increas- 
ing number of monomers. On the other, the concomitant increase of the lattice constant 
brings the neighbors further apart and therefore it slows down the growth of the repul- 
sions from the neighboring dendrimers, which have the tendency of shrinking the molecules. 
The insensitivity of the equilibrium value to the density can be understood in terms of an 
equality of the osmotic pressure from the interior of the clusters, Ilint, and that from the 
neighboring clusters, Ilext- As both scale proportionally to A^occ and the lattice constant is 
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FIG. S2. The radius of gyration of (a) bcc and (b) fee crystals as a function of crystal occupation 
number for various crystal densities as indicated in the legends. The open symbols denote the 
value of Rg{p;Nocc) at the optimal occupation number, obtained by thermodynamic integration 
as described in the main text, and they are color- and symbol coded as the filled symbols. The 
horizontal line indicates the value Rg = 3.59a of the gyration radius of the dendrimers at infinite 
dilution. 

rather insensitive to the density, this quantity cancels out from both sides of the equation, 
bringing about the aforementioned insensitivity of Rg on p. This rigidity of the dendrimers 
with respect to strong deformations lies at stark contrast to the behavior of linear or ring- 
shaped polymer chains, which shrink above their overlap density p* according to the law 

The fact that Rg{p) turns out to be also almost identical to its zero-density value does 
not imply, however, that the conformation of the dendrimer is invariant with density. In 
Fig. [S3t the monomer density profiles with respect to the dendrimer's center or mass in 



cluster crystals are shown. In particular, in Fig. S3 a) the profiles of the core monomers 
((7 = and g = 1) are compared for several occupations Nocc at fixed density pR^g = 0.282, 
including the profiles for the equilibrium occupancy, A^"^^^ = 8; as a reference, the profile of 
an isolated dendrimer is also plotted. The same comparison is shown for the shell monomers 



{g = 2) in Fig. S3 b). We find that for all finite-density profiles, and in contrast to the case 
of the isolated dendrimer, the back-folding of shell monomers to the center of the dendrimer 
is suppressed due to the presence of core monomers of other dendrimers which occupy the 
available space. Further, increasing the occupancy naturally leads to an increasing spatial 
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extension of the core- monomer aggregate as well as an increased lattice constant. Thus, the 
shell monomer distribution broadens for increasing occupancy, and the maximum probability 
distance moves slightly towards the exterior of dendrimers. It can be seen that the profiles of 
the equilibrium system (A^^^c — 8) and of an isolated dendrimer differ significantly from one 
another, while their radii of gyration were found to almost coincide. This equality originates 
in two opposing effects, namely a reduced probability of finding shell monomers within the 
dendrimer center at finite densities, paired with an increased probability of finding core 
monomers at the center. 




FIG. S3. Comparison of the monomer density profiles of dendrimers in the cluster fcc-solid for fixed 
density p and varying occupancy A'^occ, as indicated in the legend. As a reference, the monomer 
profiles of isolated amphiphilic dendrimers {pB?g = 0) are also shown, (a) Core profiles; (b) Shell 
profiles. 



[1] B. M. Mladek, G. Kahl, and C. N. Likos, Phys. Rev. Lett. 100, 028301 (2008). 
[2] D. A. Lenz, B. M. Mladek, C. N. Likos, G. Kahl, and R. Blaak, J. Phys. Chem. B 115, 7218 
(2011). 

[3] P. Welch and M. Muthukumar, Macromolecules 31, 5892 (1998). 
[4] B. M. Mladek and D. Frenkel, Soft Matter 7, 1450 (2011). 



16 



